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Abstract. We study the temperature dependence of the adjoint Polyakov loop and 
its imphcation for the momentum spectrum of gluons in the mean-field approximation. 
This allows us to calculate the contribution of the thermal (transverse) gluons to the 
thermodynamic pressure. As an application, we evaluate the rates for the strange 
quark pair-production processes qq —^ ss and gg — > ss as functions of temperature 
including thermal effects on quark deconfinement and chiral symmetry breaking. 
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1. Introduction 

The increased abundance of particles containing strange quarks in the spectrum of 
emitted hadrons, especially hyperons, was proposed by Hagedorn and Rafelski [Ij as 
a signal for the formation of quark matter in relativistic heavy ion collisions. Soon 
afterwards, the enhanced pair production of strange quarks required for the saturation 
of strange quark phase space was predicted to occur as an effect of quark and gluon 
deconfinement [21 |3l HI |5]. The predicted enhancement has been observed in many 
experiments (see, in particular: |6l [7]). The parameter 7^, describing the degree of 
saturation of the phase space of strange hadrons, has been determined by thermal 
chemical fits to the abundances of hadrons emitted in collisions between two ''Au nuclei 
at center-of-mass energies of 200 GeV/nucleon at the Relativistic Heavy Ion Collider 
(RHIC). Values for this parameter obtained in such fits range from 7^ = 1.03 ± 0.04 [8] 

to 7, = 2.00 ±0.02 ig. 

It is important to understand how the dynamics of deconfinement and chiral 
symmetry affects the prediction of production of the strange quarks and therefore its 
enhancement. The recently developed three-flavor PNJL model [TOlfTT] has made such a 
study possible. In the present work we explore the effect of the temperature dependence 
of the Polyakov loops and chiral condensates on the strange quark production. Our 
article is structured as follows. After stating the basic equations of the PNJL model, 
including the thermal quark and antiquark distribution functions, we obtain an explicit 
fit for the temperature dependence of the effective action of the Polyakov loop in mean- 
field theory. We confirm that this action satisfies the scaling of the thermal averages of 
the Polyakov loop in different color representations by their Casimir operator. We then 
calculate the temperature dependence of the thermal average of the adjoint Polyakov 
loop, the thermal distribution function of transverse gluons and the contribution of 
transverse gluons to the thermodynamic potential. 

Finally, we calculate the temperature dependence of the pair-production rate 
of strange quarks using the Polyakov loop-suppressed quark and gluon distribution 
functions. With the aid of the PNJL model, we identify the temperature where the 
gluonic contribution to the production rate becomes dominant. The cross-over of the 
contributions from light quarks and gluons is a novel phenomenon which does not exist 
in the traditional approach to strange quark-pair production based on perturbation 
theory [2]. 

2. The PNJL model 

The phase transformations of QCD matter due to deconfinement and chiral symmetry 
restoration have been combined in one theoretical framework, which is the Nambu-Jona- 
Lasinio model with the Polyakov loop (PNJL model) [I2l [13]. The Polyakov loop in 
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color-SU(3) representation r is defined as 

»1/T 



V exp ig 



(ItAa 



x,r 



(1) 



where V denotes that the exponential is path-ordered, T denotes the temperature 
and (x, r) is the temporal component of the SU(3) gauge field in representation 
r. In particular, L3 and Lg denote the Polyakov loops in the fundamental and adjoint 
representations, respectively. The traces of the Polyakov loops are defined as 

£3 = iV-^ tr^ L3, £3 = iV-i tr^ 4, (2) 
4 = (iVc -l)"'trAL8, (3) 

where tr^T^ and tr^ denote the color traces in the fundamental and adjoint representation, 
respectively. Note that the dependence of L^, Lj, Cr and Cr on the spatial coordinate x 
is suppressed in (II])-(l3]). With the above definitions, the Lagrangian of the three flavor 
PNJL model is [llj 



where D. 



^ a=0 

+ go [det (1 - 75) + h.c. , 
— g6fj_iA4^ is the gauge-covariant derivative, W(^3, 



(4) 

?3; T) is the effective 



potential for the Polyakov loop and the three-flavor current quark mass matrix tjiq 



diag^mufiyTTidflymsfi). In the limit of isospin symmetry, m„^o 



the mean-field approximation, the chiral condensates and the thermal averages of the 
Polyakov loops are the order parameters of the phase transition [HI |13]. Note that, in 
the mean-field approximation, the thermal average of the Polyakov loop is independent 
of the spatial coordinate x. 

To encode the features of the temperature dependence of the effective potential 
U{i3,i3;T), the nonperturbative contribution to the gluon thermodynamic potential 
per unit volume is assumed to be of the following form 



NP 



bT 



+ In 



54 exp 

1- 6(^3) (4) 



-f j (4) (4) 



(4)(4))' + 4((£3)^ + (4)^ 



(5) 



where (£3) and (£3) are the thermal averages of the Polyakov loops in (j2]). Standard 
values for the parameters in (jl]) and (jSD are a = 0.664 GeV, b = 0.03A'^, gs = 3.67A"^, 
= — 9.29A~^ and A = 0.6314 GeV QjJ. Moreover, it is straightforward to obtain the 
quark partition function from the quark Lagrangian, yielding the quark grand canonical 
thermodynamic potential per unit volume [TTj : 



- 2T 
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6 E / 7^^/(k)^ (A^ - IkD , (6) 



f=u,d. 
,2^1/2 



(27r)^ 

where £"/(!<;) = (|k|^ + m^)^/^ are the single-particle energies, and ruf and /i/ are the 
constituent mass and chemical potential of the quark with flavor /. The constituent 
quark masses are related to the current masses and the chiral condensates by 

mq = mq^o-2gs{qq) -2gD{qq){ss), {q = u,d) (7) 
rus = rrisfi - 2gs{ss) - 2gD{qqf, (8) 

where (uu) = (dd) and q = ^d,o- Furthermore, in the mean-field approximation, the 
quark condensate contribution to the thermodynamic potential per unit volume is [TT] 

f^cond = 9s (2(gg)' + {ssf) + 4gD{qq)\ss). (9) 

In the present article, the parameters are chosen as follows: mg_o = 0.0055 GeV, 
TJisfi = 0.1357 GeV. The vacuum quark condensates are {qq)o = (—0.246 GeV)'^ and 
(ss)o = (-0.267 GeV)^ 

In the mean-field approximation, one easily derives the temperature dependence of 
the chiral condensates (qq) and (ss) as well as the thermal averages of the Polyakov 
loops in the fundamental representation, i.e. (£3) and (£3) [TT]. Furthermore, it is 
straightforward to derive from ([6]) the color-averaged distribution functions for q and q 

m- 

{i,)X_^ + 2{l,)Xl + X^, 
l + 3(4)A+ + 3(l3)A^ + A3 

{h)X. + 2{is)Xl + X' 
l + 3(l3)A_ + 3(£3)A2 

where we introduced the abbreviations X± = exp[— ((|kp -|- m^)^/^ ip ji)/T] and chose 
an isospin-independent chemical potential for the light quarks: fi = fiu = fJ-d- 
Equations flTU]) and fITT]) are obtained by assuming the Weiss mean-field approximation 
^llj . In Appendix A we discuss several Polyakov- loop averaging procedures for 
evaluating the quark thermodynamic potential and thereby comparing the quark 
distribution functions derived from these different averaging procedures. In the 
following, we will consider the choice fi = 0.1 GeV. These distribution functions are 
to be used in the study of the strange quark pair-production rates due to the process 
qq ss. 



^"^^ l + 3(4)A+ + 3(l3)AUAr ^ ^ 

f,,,,.J§^±m^^, (n, 



3. Adjoint Polyakov loop and gluon distribution function 

In order to describe the process gg ss, we also need to calculate the adjoint Polyakov 
loop and study how it affects the gluon distribution function. Starting from the Yang- 
Mills Lagrangian Cg = —{1 / A)F'^^F°''^^ , it straightforward to write down the gluon 
partition function for the transverse gluons, and obtaining the perturbative contribution 
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to the gluon thermodynamic potential per unit volume [T5| [TB] . 

+ 8/^|k|^(A2-|kp). (12) 



. (27r)3 

Because v44(x, r) is assumed to be independent of the spatial coordinate x in the mean- 
field approximation, the Polyakov loop in the fundamental representation can be gauge 
rotated to diagonal form, L3 = diag(e*"'^^, e*"^^, e"'^^), with 03 = — (0i + 02)- Therefore, 

4 = (4)* = ^[exp(i0i) + exp(i02) + exp(-i(0i + 02))]- (13) 
In the same gauge, Lg can be expressed in terms of the eigenvalues of L3 [15j : 

Lg = diag (1,1, e'^''' , e^"^" , e'^^'' , e'^''^ , e'^^' , e'^^^ ) , (14) 
where 0^^ = 0^ — 0^. Thus, 

4 = ^tr^Ls = ^ I 1 + Zlcos0jfc j • (15) 

\ j<k / 

Inserting f[T^ into f[T^ and after some algebraic transformations, the color-averaged 
gluon distribution function is obtained as 

/. (^) = o E (t^A L^) exp (-n |k| /T), (16) 



n=l 



where 



tr^ 



1 + ^ cos ncpjk 
j<k 



[17) 



The remaining task is to evaluate (tr^ Lg). In order to do so, we need to specify the 
full distribution of eigenvalues of the Polyakov loop, i. e. the distribution of phases 0j. 
The thermal average of any function of the eigenvalues of the Polyakov loop, -F(0i, 02), 
is defined as 

A))) = /o^^ ^'^1 /o^^ d(j)2H{(j)i, 02)W^(01, 02; T)F{(j)i, 02) 
^ ^ /o2'^rf0i/o2'^rf02if(0i, 02)1^(01, 02; T) ' ^ ^ 

where if(0i,02) is the SU(3) Haar measure, the weight function Vr(0i,02;T) denotes 
the distribution of eigenvalues of the Polyakov loop. The Haar measure for the SU(3) 
symmetry group is given by if(0i, 02) = nj<fc sin^(0ifc/2). Since 03 = -(0i + 02), the 
integration in (JTSl) only goes over 0i and 02. A suitable choice of the weight function is 
a crucial step in evaluating (-F(0i, 02)). We follow Gocksch and Ogilvie [17] and Gupta 
et al. [H] in choosing a weight function of the form 

1^(01, 02; T) = exp (6 rf/?3 (4) Re(4)) , (19) 

where c? = 3 and l3z{T) is a fit parameter depending on temperature. The particular 
form (fT9l) is suggested by the strong coupling expansion of the gauge theory. We note 
that exp[— ^(£3, ^3; T)/(5T)] and if(0i, 02)IV(0i, 02) have a corresponding structure 
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Figure 1. The fit parameter as a function of temperature. 

expressed in terms of the eigenvalues of the fundamental Polyakov loops, as can be seen 
from assuming ^(£3, £3; T) to be in the form of (5) with (£3) and (£3) replaced by £3 and 
£3 respectively. Starting from (|T9ll and inserting the known values of (£3) into ( ITSi) . the 
temperature dependence of jS^ can be solved numerically, as shown in figure [TJ We now 
have obtained an explicit expression for Vr(0i,02) at each temperature. 

We note that the temperature corresponding to the minimum of (3^ in figure [1] 
coincides with the critical temperature of the deconfinement phase transition in the 
mean-field approximation. In their investigation of the gluonic contribution to the 
thermodynamic potential of the PNJL model, Megias et al. [16] used a weight 
function of similar form as ( fT9i) . but did not make the mean- field approximation. They 
determined the parameter /?3 from the empirical relation between string tension and 
the deconfinement temperature of the pure gauge theory. Here, we have determined 
by imposing a self-consistency condition on the expectation value of the fundamental 
Polyakov loop. 

With the weight function (|T9i) . equation ( |T8l) allows us to evaluate the thermal 
average of the adjoint Polyakov loop and thus the gluon distribution function fg{k). 
The temperature dependence of the quark condensates and the thermal averages of 
the Polyakov loop (^3), (£3) were first calculated in [11], as well as the temperature 
dependence of (^s) are shown in figure O Furthermore, the temperature dependence of 
the constituent masses of quarks are obtained from ([7]) and ([8]), as shown in figure [31 

To verify the validity of (£3) obtained by this procedure, we check its consistency 
with the Casimir scaling of the thermal averages of the Polyakov loop observed in lattice 
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Figure 2. The temperature dependence of the order parameters, £3, £3, is, 
{qq)/{qq)o and {ss)/{ss)o. The values of the vacuum quark condensates are {qq)o = 
(-0.246 GeV)^ and (ss)o = (-0.267 GeV)^ 



QCD [18], [19]. Casimir scaling refers to a relation, valid for all temperatures, between 
the thermal averages of the Polyakov loop in different representations r of color-SU(3) 
of the form 

(4) = (4)'% (20) 

where dr = C2(r)/C2(3) and C2(r) denotes the eigenvalue of the quadratic Casimir 
operator in representation r. For the adjoint representation, ds = C2(8)/C2(3) = 9/4 
and thus {is) = {hY^'^- Figure H] shows our results for together with {i^)^^^, as 
function of (£3). As can be seen, the values of {is) and (£3) show a good agreement with 
Casimir scaling (120|) for temperatures T ranging from 0.12 GeV to 0.30 GeV. 

By numerically evaluating (tr^ Lg ) for each power n and as a function of 
temperature T, the gluon distribution function can now be obtained from f|T6l) . The 
ratios fq{k)/f-pT){k) and fg {k) //be (k) are plotted in figure [5], where /fd and /be are 
the free Fermi-Dirac and Bose-Einstein distributions, respectively, with the constituent 
quark mass niq in ((Tj) and /i = 0.1 GeV. When T is near T^, gluon is more strongly 
suppressed than quarks in spite of the additional effect of chiral symmetry breaking on 
the constituent quark mass irtq. Furthermore, we evaluate the quark, antiquark, net 
quark and gluon number densities by 

ng = nu + nd = ANcJ-^^fgik), (21) 
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Figure 3. The temperature dependence of the constituent quark masses, rriq and 
in ([7]) and ([S]) respectively. 



nq = nu + = 4iVc / J^^hi^). (22) 



^ im) - f,{k)] , (23) 



n, = 2(iV,^-l) J j^JAk), (24) 

where the momentum integration is taken without imposing any cutoff. The number of 
ffavors in (!2T!) - (l23l) is 2 because we are evaluating the quark, antiquark, or net quark 
number densities for flavor u and d. Figure [6t^a) shows the temperature dependence 
of the (scaled) quark, antiquark, net quark and gluon number densities, Ug/T^, rig/T^, 
Hq-g/T^ and Ug/T^, respectively. In figure [6](a), ng_g/T^ possesses similar features 
as that in the two- flavor PNJL model [13]. Moreover, Ug/T^ is non- vanishing at low 
temperature due to the existence of the color-singlet gluon states. On the other hand, 
assuming zero quark chemical potentials for all three flavors, yU = /i^ = 0, we evaluate 
Ug/T^, Hg/T^ and the (scaled) strange quark number density rig/T^ as functions of the 
temperature, as shown in figure [6]J^b). We note that at T = 1.5 the quark densities 
have reached almost 90% of their asymptotic values; whereas the gluon density is still 
less than 2/3 of its asymptotic value (n^/T^ ^ 1.95). The quark and gluon number 
densities reach 99% and 92% of their asymptotic values respectively at T = 3.5Tc. 

With fg{k), we can now calculate the contribution of the thermal (transverse) 
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Figure 4. The relations of (Is) to {£3) and (^3}^^'' to (is). The temperature T ranges 
from 0.12 GeV to 0.30 GeV. These two curves indicates that the values of {£s) and 
(^3} is consistent with Casimir scaling, {ig) — {£3}^^^. 
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Figure 5. Ratios of distribution functions, (a) fq{k)/ fFuik), with u and d quark 
chemical potential ^ = 0.1 GcV, and (b) f g{k) / f-Q^{k) , at temperatures, T — 0.3, 0.2 
and 0.18 GeV. 



^luons to the thermodynamic pressure: 



Pg 



(25) 
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Figure 6. (a) The temperature dependence of the (scaled) quark, antiquark, net quark 
and gluon number densities, Uq/T^, riq/T^^ Uq-q/T^ and Ug/T^, respectively. The u 
and d quark chemical potential = 0.1 GeV. (b) Ug/T^, riq/T'^ and the (scaled) strange 
quark number density Ug/T^ as functions of the temperature, under the assumption 
of zero quark chemical potentials for all three flavors, = = 0. In (b), the three 
dashed curves represent asymptotic values of Ug/T^, '^ql^^^ s-nd f^s/T^ respectively, 
for {t^) = (£3) = (^g) = 1 and the current quark masses. 



The (scaled) pressure PglT^ is plotted as a function of the temperature in figure [71 
Because the thermal average of the adjoint Polyakov loop {^^) 1 for T ^ T^., 
Pg/T^ approaches the asymptotic behavior predicted by the Stefan-Boltzmann value, 
Pg/T^ 167r^/90 ~ 1.75 when T ^ T^. Figure [7] confirms that the gluon pressure 
reaches 95% of the Stefan-Boltzmann value at T = 700 MeV. 



4. Strange quark pair-production rate 

With the expressions for fg{k), fq{k) and fg{k), the strange quark pair-production rate 
for both qq ss and gg — » ss can be derived in the PNJL model. The strange quark 
pair-production rate per unit volume is given by 

where 



Aq = - dss\ 1 ^ 6(s-{ki + k2) ) aqg^ss (s) 



2 J4m1 



1 /-oo . 2\ _ 

(27r)3|ki| i (27r)3 jkal V2 



gg^ss \S) 
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Figure 7. The temperature dependence of the (scaled) thermodynamic pressure, 
Pg/T*, contributed by the thermal (transverse) gluons. 



where Eg(k) = {\kf + nilf/^. The cross sections are exphcitly given by 



27^ 



s +2s [niq + 1 + Wm^m', 
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1 + M + ^jtanh 
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Setting ki = |ki| and k2 
9 



\« 8s J \ s J 
|k2|, we can simphfy f l27|) and 



ds s \ 1 



Arrfi 
a 



qq- 



as follows: 

kik2 



X 9 
A 



2 (klk2 + Eg{kl)Eg{k2) + 







fq{kl)fq{k2 



^q — ~^ / C^-S S Cqg^ss 



dkidk2 9 {Akik2 - s) fg (ki) fg (^2^ 



(29) 



(30) 



(31) 
(32) 



By substituting equations ([TO]), ([II]), ([16]), ([29]) and ([30]) into equations ([3l]) and ([32]), the 
numerical values of the production rates are obtained as functions of the temperature. 

Figure [8] shows the temperature dependence of the strange quark pair-production 
rates in the PNJL model, compared with those obtained for free quarks. One notices a 
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number of qualitative differences between the rates calculated in the the PNJL model 
and those calculated in free perturbation theory. First, the rates are suppressed for all 
values of the temperature. This is, in part, due to the suppression of the thermal quark- 
and gluon excitations by the Polyakov loop and, in another part, due to the fact that 
the effective strange quark mass remains larger than the current quark mass even at 
temperatures moderately above T^, as shown in figure O We also note that the curves for 
the Polyakov loop-suppressed gluon induced production rate Ag drops below the quark 
induced production rate Ag below T ^ 240 MeV, reflecting the stronger suppression of 
gluons at low temperature caused by the Casimir scaling of the thermal average of the 
adjoint Polyakov loop. 

5. Conclusions 

We have studied the effects of deconfinement and chiral symmetry breaking on the rates 
of strange quark pair production in the framework of the PNJL model. As proposed in 
[2], H], the strange quark pair-production rate is enhanced in the deconfined phase for 
the free quarks and the production rate for gg — > ss is dominant at all temperatures. In 
the PNJL model, the enhanced production of strange quarks is also obtained, but the 
production rates for qq ss and gg — > ss cross over at ~ 240MeV. The production 
rate for qq ^ ss is dominant when T < Tr, while that for gg ss is dominant when 
T > Tr. Besides, when T < Tc, the production rates for qq — > ss and gg ss are 
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both very small in the PNJL model because quark and gluon quasiparticles are strongly 
suppressed below Tc. In this temperature region, strangeness production is dominated 
by hadronic reactions, which were investigated by Rehberg et al. in the NJL model |20j . 

In figure [HI the coupling scales out when we compare the quark and gluon 
contributions to the production rate. We note that, approaching Tc, one needs to take 
into account the interactions originated from the appearance of collective modes due 
to the onset of the spontaneous breaking of chiral symmetry, rendering our treatment 
incomplete in the transition region. However, our goal was to study at what temperature 
above the gluonic contribution to the production rate becomes dominant. We found 
that this temperature is around 240 MeV within the framework of the PNJL model. 
Because this threshold is well beyond the temperature range in which the chiral phase 
transition occurs, as can be seen from figure [21 our neglect of the contribution from 
collective (hadronic) modes appears justified. 

A by-product of our investigation is the demonstration that the thermal average in 
( ITSl) satisfies Casimir scaling (1201) of the fundmental and adjoint Polyakov loops. This 
gives confidence that the weight function f[T9|) can be used to obtain the temperature 
dependence of any quantity that involves the eigenvalues of the Polyakov loop. For 
example, the temperature dependence of fg{k) requires the evaluation of the averages 
(tr^Lg). In turn, fg{k) makes it possible to compute the contribution of the thermal 
(transverse) gluons to the pressure as a function of the temperature. This makes it 
unnecessary to include this contribution explicitly in the effective potential for the 
Polyakov loop, as sometimes done in the literatures [2T]. 

In our study, the phase transformations of QCD, including deconfinement and chiral 
symmetry breaking, are incorporated into the evaluation of the thermal strange quark 
pair-production rate. Using the same techniques, the effects of the Polyakov loop on 
other signatures of quark-gluon plasma can be explored in the future. 
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Appendix A. Averaging procedures 

In this section, we discuss several Polyakov-loop averaging procedures for the grand 
canonical thermodynamic potential in the quark sector and study their implications for 
the quark and antiquark distribution functions. This study can be easily extended to 
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the Polyakov-loop averaging procedures in the gluon thermodynamic potential, which 
are not exphcitly formulated here. The quark grand canonical thermodynamic potential 
per unit volume is defined in terms of the quark grand partition function, 

Q^ = -lln{Z), (A.l) 

where V denotes the volume of the system. In the mean-field approximation, the quark 
grand partition function is associated with a set of quantum numbers, a = {k, s, f, c, ±}, 
where k, s, f, c and ± denotes momentum, spin, fiavor, color and particle/antiparticle 
quantum number respectively. The average in flA.ll) is taken over the eigenvalues of the 
Polyakov loop, as shown explicitly in ( ITSj) . Instead of the full average used in ( lA.ll) . an 
approximate averaging method used more frequently in the PNJL model is the quenched 
average: 

fi,^ -^(InZ) = -^(IndetZ,), (A.2) 

where denotes the single-particle partition function for each quantum number and 
the determinant runs over all quantum numbers a. Various further approximations 
can be applied to ( lA.2p . which entail distinct averaging procedures. In the following 
texts, we discuss the differences among several Polyakov-loop averaging procedures for 
the quark thermodynamic potential and clarify their effects on the quark and antiquark 
distribution functions. 

Define the following subsets of a: a = {k, s, /, ±}, ai = {k, s, /} and 0:2 = 
Starting from ( 1A.2I) . we have the following (approximate) averaging procedures: 

^1- -^Eln(detZ^A, (A.3) 
-^E(liidetZ^,A, (A.4) 

-^Eln(ndet^<,„<,„c), (A.5) 

where detc denotes the color determinant. Equation (lA.SP is the Weiss mean-field 
approximation, which is frequently used in the literatures of the PNJL model [HI [T3]. 
We note that flA.SP and flA.4|) take the Polyakov-loop average for quarks and antiquarks 
separately. This implies that, in the limit {tip L^) —>■ 0, only states with baryon quantum 
(quark-triplet) numbers contribute, but not states with meson quantum numbers (quark- 
antiquark pairs). We further note that (1A.3I) and ( ]A.5l) replace the quenched average 
over the Polyakov loop configuration by the unquenched average. We compare Qq with 
flq in Appendix A.l and Qq with Clq in Appendix A.2 respectively, by deriving the 
quark distribution functions from equations ( ]A.3P -( lA75l) . 
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Appendix A.l. Validity of the Weiss mean- field approximation 



By the averaging procedure in 
volume in ([6]) is simplified to be 



3]), the quark thermodynamic potential per unit 



2T y: 

f=U,d,! 



d^k 



{in 



(27r)3 

+ ln [l + 3 (4) A- + 3(£3)A'_ + A^ 



l + 3(£3)A+ + 3(l3)A^ + A^ 
} 



-6 E 

f=u,d,s 



(0W''(A'-Ik| 



(A.6) 



where A± = exp[— ((|k| + nig) ' =l= lj)/T]. The quark and antiquark distribution 
functions, (fTOl) and (fTTj) . are easily obtained from (]A.6p . On the other hand, if we start 
from ^ and use the averaging procedure defined in (lA.4p . the quark and antiquark 
distribution functions are alternatively obtained: 

/ 4A+ + 2£3A^ + XI \ 



1 + 34A+ + ^h>A + A^ 
4A- + 2£3A^_ + A^ 
34A2 + A3 



(A.7) 



(A. 



\ 1 + 3£3A_ -r 0^3/ 

where £3 and £3 are expressed in (fT3il . By the definition of the thermal average (JT8l) and 
the weight function f|T9|) . we can evaluate (IA.7I) and flA.SP explicitly. Without losing 
generality, in this section we assume a vanishing u and d quark chemical potential, i.e. 

= 0, which implies the simplification (£3) = (£3). Figure RT] (a) shows the comparison 
of the ratio fq{k)/f-p-Q{k) obtained from (ITOl) with the ratio JqiJ^) I f-pnik) from (]A.7p . The 
figure shows that the two ratios agree well for all temperatures, especially in the high 
momentum region, which is most relevant for the thermal strange quark pair-production 
rate. 

The gluon distribution function in (fT6|) is derived from the gluon thermodynamic 
potential obtained by the averaging procedure analogous to (lA.4p . On the other hand, 
if using the averaging procedure analogous to ( 1A.3I) instead, we obtain the gluon 
distribution function under the Weiss mean-field approximation, namely 

|(tr^L8)exp(-|k|/T) 



l-i(tr^L8)exp(-|k|/T)^ 



(A.9) 



assuming 

{iiALD ^ % {^{ii aL,)^ . (A. 10) 

Figure IXTj fb) shows the comparison of the ratio /c,(/c)//be(^) obtained from (lA.Op with 
that from f|T6|) . They agree well in the high momentum region but deviate in the low 
momentum region. 



Appendix A. 2. Alternative quark and antiquark distribution functions 

By the Polyakov-loop averaging procedure (lA.Sp . we evaluate alternative quark and 
antiquark distribution functions, which, unsurprisingly, incorporate the probabilities of 
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Figure Al. Validity of the Weiss approximation for, (a) fq{k)/ fpriik), assuming a 
zero quark chemical potential, and (b) fg{k)/fBE{k)- The solid curves denote the 
ratios obtained by, (a) /q(fc)//FD(fc), where fq{k) is defined in (|A.7|) with ^ = 0, 
and (b) /g(fc)//BE(fc), where fg{k) is defined in (fT6)l . The dashed curves denote the 
ratios obtained by, (a) f q{k) / f-p-oik) , where fq{k) defined in (fTO)l with ji — Q and 
(£3) = (£3), and (b) fg{k)/ fsEik), where fg{k) defined in (|A.9p . The curves are 
obtained at temperatures, T = 0.3, 0.25, 0.2, 0.18 and 0.15 GeV (top to bottom). 



color-singlet quark-antiquark states. The quark distribution function derived from the 
averaging procedure (lA.SP is 



3e^^ + 3{tTF L3)e t + 2{tiF Ll)e t 



+2(trj. L3 tip Ll)e^^ + 3(trF Ll)e ? 

-(4fiq(fc)+2M) ^ , > , , -(3Eg(fc)-M) 

+ (tr^ L3)e T + 2(2(tr^ L3) + (tr,.(4)2))e ^ 

+ {2{trFLl) + {tip Ll))e t + (tri^Latr^. Ll)e^^ 
+3e T + 2(tri7 Ll)e t + (tr^ 1^)6 t 



X 



~6Eq{k) -(5Eq(fc)-M) , X, -(5gg(fc)+M) 

1 + + (tr^ L3)e t + (tr^ Ll)e 



T 



, ~4fig(fc) , -(4Eq(fc)-2M) 

-|-(tri7 L3 triT Lg/e t + {tip L^'i)'^ ^ 

+ (tr,.L3)e ^ + (2(tr^L3) + (tr^(4)2))e t 

+ (2(tr^4) + (tr^L2))e t 

, -2Eq{k) -3(Eg(fc)-M) -3(gg(fc)+M) 

+ (trirL3trirL3)e t +e r +e t 

, -2(Eg(fc)-M) , ^ -2(Eq(fc)+M) 

+ (tr^LT)e T + (tr^L3)e t 

-(Eg(fc)-M) , . -(Eq(fc)+M) l 

+ (tr^L3)e ^ +(trpLT)e ^ , (A.ll) 
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k (GeV) 



Figure A2. The comparison of tlie ratio /g(fc)//FD(fc) obtained from (jA.lip . in solid 
curves, with fq{k) / fpuik) obtained from (jlOp . in dashed curves, assuming a quark 
chemical potential /i = 0. The curves are obtained at temperatures, T — 0.3, 0.25, 0.2, 
0.18 and 0.15 GeV (top to bottom). 



where Eq{k) = (fc^ + mg)^/^ and {tip Lstrp lI) = (tiyiLg) + 1. Moreover, the antiquark 
distribution function fq{k) can be obtained from (lA.llI) by interchanging L3 and L3 
and replacing /x by — /i. Equation (lA.lip contains a sum over the probabihties of all 
states of Ni quarks and N2 antiquarks, where Ni,N2 = {0,1,2,3}. As noted before, 
(lA.llI) incorporates the contribution of color-singlet quark-antiquark states, which are 
not contained in the expression ( fTOl) . Figure IA2I shows the comparison of the ratio 
fq{k) / fY-oik) obtained from flA.lip with fq{k)/f-pT>{k) obtained from flTUl) . again for 
/i = 0. The figure shows that the two ratios agree well for all temperatures, especially 
in the high momentum region. We also note that, when T is near T^oi T < T^, fq{k) in 
(lA.llI) is slightly larger than fq{k) of ffTOl) in the low momentum region. This difference 
can be traced back to the contribution of the color-singlet quark-antiquark states. 

In conclusion, because both fq{k) and fq{k) are in good numerical agreement with 
fq{k), we are justified to use (fTOj) and (ITT!) in the evaluation of the strange quark pair- 
production rate, as done in the main part of this article. 



References 



[1] R. Hagedorn and J. Rafclski, Phys. Lett. B 97, 136 (1980). 

[2] J. Rafclski and B. Miiller, Phys. Rev. Lett. 48, 1066 (1982) [Erratum-ibid. 56, 2334 (1986)]. 
[3] T. S. Biro, B. Lukacs, J. Zimanyi and H. W. Barz, Nucl. Phys. A 386, 617 (1982). 



18 



[4] P. Koch, B. Miiller and J. Rafelski, Phys. Kept. 142, 167 (1986). 

[5] T. Matsui, B. Svetitsky and L. D. McLerran, Phys. Rev. D 34, 783 (1986) [Erratum-ibid. D 37, 
844 (1988)]. 

[6] E. Andersen et al. [WA97 Collaboration], Phys. Lett. B 449, 401 (1999). 

[7] B. I. Abelev et al. [STAR Collaboration], Phys. Rev. C 77, 044908 (2008). 

[8] M. Kaneta and N. Xu, arXiv:nucl-th/04 05068, 

[9] J. Letessier and J. Rafelski, Eur. Phys. J. A 35, 221 (2008). 

[10] M. Ciminale, R. Gatto, N. D. Ippolito, G. NarduUi and M. Ruggieri, Phys. Rev. D 77, 054023 
(2008). 

[11] K. Fukushima, Phys. Rev. D 77, 114028 (2008) [Erratum-ibid. D 78, 039902 (2008)]. 

[12] K. Fukushima, Phys. Lett. B 591, 277 (2004). 

[13] C. Ratti, M. A. Thaler and W. Weise, Phys. Rev. D 73, 014019 (2006). 

[14] H. Hansen, W. M. Alberico, A. Beraudo, A. Mohnari, M. Nardi and C. Ratti, Phys. Rev. D 75, 
065004 (2007). 

[15] P. N. Meisinger, T. R. Miller and M. C. Ogilvie, Phys. Rev. D 65, 034009 (2002). 

[16] E. Megi'as, E. Ruiz Arriola and L. L. Salcedo, Phys. Rev. D 74, 065005 (2006). 

[17] A. Gocksch and M. Ogilvie, Phys. Rev. D 31, 877 (1985). 

[18] S. Gupta, K. Hiibner and O. Kaczmarek, Phys. Rev. D 77, 034503 (2008). 

[19] S. Gupta, K. Hiibner and O. Kaczmarek, Nucl. Phys. A 785, 278 (2007). 

[20] P. Rehberg, S. P. Klevansky and J. Hiifner, Phys. Rev. C 53, 410 (1996). 

[21] S. Roessner, C. Ratti and W. Weise, Phys. Rev. D 75, 034007 (2007). 



